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'This  report  provides  the  results  of  a study  of  attenuation  and  depolarization  of 
millimeter  waves  along  a satel 1 i te-to-earth  communication  Jink.  Specifically, 
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2)  Development  of  physical  models  of  rain  attenuation  and  depolarization.  ^ 
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FOREWARD 


The  goals  of  this  research  project,  as  detailed  in  the  proposal. 


were  to  provide  additional  support  for  the  study  of  attenuation  and  depolari- 


zation of  millimeter  waves  along  a satellite-to-earth  communication  link. 


Specifically,  this  support  contributed  to  our  research  efforts  in 


1)  Detailed  analysis  of  attenuation  and  depolarization  data  from  the  CTS 
(Communications  Technology  Satellite)  satellite  downlink  beacon  at  a 


frequency  of  11.7  GHz,  and  the  19.04  and  28.56  GHz  downlink  beacons  on 


the  COMSTAR  satellites,  and  2)  Development  of  physical  models  of  rain 


attenuation  and  depolarization.  These  two  activities  are,  of  course,  very 


closely  related.  It  is  necessary  to  have  a reliable  data  base  for  the 


development  of  a theoretical  model.  On  the  other  hand,  data  collected 


for  a specific  location  frequency,  polarization,  and  elevation  angle 


cannot  be  extrapolated  to  other  locations,  frequencies,  polarizations,  and 


elevation  angles  with  full  confidence.  Our  approach  is  to  build 


a model  which  accurately  predicts  propagation  effects  as  verified  by 


comparison  to  data  for  several  locations,  frequencies,  polarizations  and 


elevation  angles,  and  then  use  the  model  to  predict  attenuation  and 


depolarization  for  situations  where  no  data  exists  but  applications  do. 


A reliable  model  eliminates  the  time  and  expense  of  experimental  data 


collection  in  advance  of  an  operational  system.  It  is  the  aim  of  this 


research  effort  to  construct  such  a data -confirmed  theoretical  model. 


Since  this  research  has  received  multiple  agency  support,  documentation 


of  the  results  Is  found  in  several  sources.  Very  detailed  reports  on  the 
data  collected  are  found  in  [l]-[3].  Information  on  the  experimental 
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hardware  Is  found  In  [4].  Details  on  the  data  reduction  procedure  used  are 
contained  In  [1], 

Support  from  the  U.S.  Army  permitted  concentration  of  effort  on  an 
existing  theoretical  analysis  effort  in  the  form  of  a Rain  Propagation 
Prediction  (RPP)  computer  program.  This  program  was  used  to  develop  a 
Synthetic  Storm  Model  (SSM).  The  synthetic  storm  model  provides  data 
concerning  the  physical  nature  of  the  rain,  such  as  the  rain  rate  distri- 
bution, rain  path  extent,  and  rain  drop  canting  angle  distribution. 

This  together  with  system  parameter  values  (frequency,  elevation  angle,  and 
antenna  polarization)  is  used  as  input  with  the  RPP  program  to  generate 
predicted  values  of  signal  fade,  isolation,  and  phase  as  a function  of 
earth  station  rain  rate.  Details  of  this  theoretical  model  development 
are  found  In  [5]. 

In  this  report  recent  results  and  essential  conclusions  of  the 
modeling  investigation  are  presented.  There  are  two  major  findings  which 
are  directly  traceable  to  the  theoretical  effort  support  by  the  Arrny: 
the  synthetic  storm  model  and  the  nonreciprocal  effects  of  rain  on 
propagation.  These  two  topics  are  discussed  in  this  report. 


! i 
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A SYNTHETIC  STORM 


MODEL  FOR  PREDICTION 


OF  MILLIMETER  WAVE 


RAIN  ATTENUATION 


A.  INTRODUCTION 


Theoretical  procedures  for  predicting  millimeter  wave  attenuation 
Introduced  by  rain  usually  require  a uniform  rain  rate  along  the 
propagation  path.  However,  nonuniform  rain  rate  conditions  frequently 
occur  on  earth-satellite  paths.  The  total  attenuation  introduced  by 
real  (nonuniform)  rain  is  of  Interest  for  communication  system  applications. 
The  effective  path  length  approach  is  a popular  method  for  predicting 
total  attenuation  on  a satellite  path,  but  we  will  show  that  the 
effective  path  length  formulation  is  Inherently  restrictive. 

Total  attenuation  can  be  determined  more  directly  by  the  synthetic 
storm  model,  which  is  introduced  here.  This  model  includes  the 
nonuniform  character  of  the  rain  in  a way  that  facilitates  prediction  for 
various  elevation  angles  and  frequencies.  Attenuation  data  from  several 
experiments  are  presented  in  support  of  the  model.  Scaling  of 
attenuation  values  from  one  frequency  to  another  is  discussed.  Also, 
the  synthetic  storm  model  together  with  the  Rice-Holmberg  rain  rate 
equation  [6]  is  used  to  predict  attenuation  statistics  from  rainfall 
accumulation  data. 
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B.  THE  EFFECTIVE  PATH  LENGTH  MODEL 


Many  of  the  theoretical  developments  in  rain  attenuation  prediction 
have  centered  on  specific  attenuation  a,  which  is  the  signal  attenuation 
caused  by  a uniform  rain-filled  space  one  kilometer  in  extent.  Specific 
attenuation  is  a function  of  both  the  rain  and  system  parameters  and 
can  be  expressed  by  the  power  relation 

bl 

a(R,f,e)  = a1  R 1 dB/km  (1) 

where  R is  the  rain  rate  in  mm/hr,  f is  the  frequency,  and  e is  the 
elevation  angle  of  the  propagation  path.  Constants  a-j  and  b-j  are  available 
from  1 to  1000  GHz  [7];  however,  in  [7]  they  were  derived  assuming  spherical 
raindrops.  A more  accurate  drop  shape  distribution  [5]  uses  a mixture  of 
oblate  spheroidal  and  spherical  raindrop  shapes.  Although  not  shown 
explicitly  in  (1),  a is  a function  of  the  fraction  of  oblate  drops  as  well 
as  the  polarization  state  of  the  incident  wave  and  the  orientation  of  the 
electric  field  relative  to  the  canting  angles  of  the  oblate  drops.  The 
parameter  dependences  of  a must  be  considered  when  developing  a total 
attenuation  model  and  particularly  when  a power  law  fit  is  employed. 

The  total  attenuation  due  to  a rain  of  physical  extent  L along  the 
propagation  path  is  given  by 

L 

A(f,e)  = / o{R(0,f,e}  d£  dB  (2) 

0 

where  «.  is  the  position  along  the  rain  portion  of  the  path.  Usually  the 
rain  rate  as  a function  of  position  is  not  known,  therefore,  (2)  is 
not  evaluated  directly.  Instead,  the  effective  path  length  model  is  used 
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frequently.  Effective  path  length  Is  defined  as 


Le^Rav,f,e^ 


km 


(3) 


where  the  average  rain  rate  is  given  by 
! L 

Raw  = r / R(0  d£  irni/hr  . 

av  l o 


(4) 


The  effective  path  length  is  obtained  from  (3)  using  experimental  values 
of  total  attenuation  A and  theoretical  values  of  specific  attenuation 
a,  and  matching  values  of  attenuation  and  rain  rate  observed  for  equal 
percentages  of  time.  Typical  effective  path  length  values  are  shown  in 
Fig.  1.  These  data  were  collecteld  at  VPI&SU  along  the  nearly  same  rain 
paths.  [2]  The  19.04  GHz  data  is  for  the  calendar  year  1978  and  is  a 
combination  of  measurements  taken  with  the  COMSTAR  D-2  and  0-3  satellite 
beacons  at  elevation  angles  of  44°  and  46°,  respectively.  The  28.56  GHz 
data  were  collected  during  April,  May,  and  June  1978  from  the  COMSTAR  D-2 
satellite  at  an  elevation  angle  of  44°.  It  is  evident  from  Fig.  1 that 
effective  path  length  is  a strong  function  of  rain  rate.  In  this  case 
the  frequency  dependence  is  rather  mild;  however,  in  general  L0  will 
depend  on  frequency;  see  [5,8].  The  frequency  dependence  arises  from 
the  implicit  assumption  contained  In  (3)  that  the  rain  is  uniform  at  a 
rate  of  Rav  over  a path  of  length  Lg.  If  this  were  true  then  effective 
path  length  and  physical  rain  extent  would  be  the  same,  Le=L,  and  Lg  would 
not  be  a function  of  frequency.  The  frequency  dependence  of  A and  a does 
not  cancel  In  (3)  because  the  drop  size  distribution  (and  thus  attenuation) 
changes  with  rain  rate  in  a nonlinear  fashion.  Thus  different  rain  rate  dis- 
tributions with  the  same  average  rain  rate  will  lead  to  different  total 
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EFFECTIVE  PATH  LENGTH  (km) 


• 19  GHz 
x 28  GHz 


40  60 

RAIN  RATE  (mm /hr) 


Effective  path  lengths  for  the  VPI&SU 
19.04  and  28.56  GHz  COMSTAR  experiments 
(B8  and  CIO). 


I 

KM 


fl  H 


attenuations,  while  specific  attenuation  remains  unchanged  because  It  Is 


based  on  a uniform  rain  rate. 


For  a particular  average  rain  rate  (obtained  by  averaging  over  a 
rain  gauge  network  or  by  a single  rain  gauge),  frequency,  and  elevation 


angle,  effective  path  length  can  be  used  to  predict  attenuation  from 


(3)  as  A=a  Lg.  However,  specific  attenuation,  and  thus  effective  path 


length,  depend  upon  the  wave  polarization  state  and  the  fraction  of 


oblate  drops.  Variations  of  up  to  5 dB  in  the  predicted  total  attenuation 


can  result  if  the  correct  system  parameters  are  not  used  when  determining 


a^  and  b-j. 


C.  THE  SYNTHETIC  STORM  MODEL 


In  this  section  a quasi -physical  model  of  real  rains  is  developed 
to  eliminate  the  need  for  effective  path  lengths.  The  model  accounts  for 
nonuniform  rain  rates  directly  and  is  based  on  a series  of  simplifying 
assumptions  which  are  justified  by  agreement  with  experimental  results. 
Based  on  the  data  of  Fig.  1,  two  assumptions  will  be  made:  1)  rain  rate 
is  uniform  for  low  rain  rates,  and  2)  as  average  rain  rate  increases, 
the  rain  rate  distribution  R(t)  becomes  increasingly  nonuniform. 

The  first  step  in  the  development  of  the  synthetic  storm  model  is 
the  establishment  of  the  effective  storm  extent.  Let  the  rain  be  of 
height  H (above  a flat  earth)  and  basal  length  B (in  a plane  containing 
the  line-of-sight  path  and  local  vertical  at  the  earth  terminal  location). 
The  storm  height  H follows  from  the  0°C  isotherm  height.  [9]  For 
simplicity  we  shall  divide  the  U.S.  into  site  latitude  classes  giving 


3.5  km 

High-latitude  (above  40°) 

4.0  km 

Mid-latitude 

(5) 

4.5  km 

Low-latitude  (below  33°) 

Confidence  in  these  height  values  can  be  gained  by  examining  attenuation 
data  from  several  experiments  at  different  locations,  elevation  angles, 
and  frequencies.  Using  attenuation  values  at  10  mm/hr,  where  the  rain  is 
assumed  to  be  uniform,  the  rain  path  length  may  be  determined  as  L=A(R=10 
mm/hr)/a,  and  a is  the  calculated  specific  attenuation  at  10  mm/hr  for  that 
elevation  angle,  frequency,  polarization,  etc.  So  for  each  experiment 
there  is  a value  of  L and  e,  as  indicated  in  Fig.  2 for  a few  of  the 


Rain  extent  model  developed  from  the  synthetic  storm  model. 


experiments  listed  in  Table  I.*  These  points  also  serve  as  a guide  In  the 
selection  of  the  storm  basal  extent.  The  value 


B = 10.5  km 


(6) 


fits  the  points  Cl  and  B1  well  as  shown  in  Fig.  2.  This  is  a more 
realistic  model  than  the  common  assumption  of  infinite  horizontal  rain 
extent.  To  summarize  the  storm  extent  model,  the  rain  path  length  along 
a propagation  path  with  elevation  angle  e is 


L = 


H CSC  e e ^ eQ 


B sec  e e i e0 


(7) 


-1  H 

where  H and  B are  given  in  (5)  and  (6)  and  eQ  is  tan  g which  is  23.2°, 
20.9°,  or  18.4°  for  low,  mid,  and  high  latitude  classes. 

The  nonuniform  nature  of  the  rain  is  modeled  by  dividing  the  rain 
path  L into  N equal  intervals,  each  of  which  contains  rain  of  uniform 
rain  rate,  as  shown  in  Fig.  3.  Cell  1 is  near  the  ground  antenna  of  an 
earth-space  link.  The  rain  rates  of  all  cells  are  related  to  the  rain 
rate  of  cell  1,  R^,  by  the  following 


R 


i 


R.|  <_  10  mm/hr 


R.j  >_  10  mm/hr 


(8) 


where  {x^}  are  real  numbers.  This  relationship  includes  the  uniform 


J 


* In  this  report  the  experiments  listed  in  Table  I will  be  referred  to 
by  the  code  numbers,  i.e.  A8,  CIO,  etc. 
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Figure  3.  Two  level  rain  rate  model  for  the  synthetic 
storm  model;  see  (10). 
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rain  rate  assumption  for  rates  under  10  mm/hr  and  a piecewise  uniform 
rain  model  for  rates  above  10  mm/hr.  The  total  rain  attenuation 
corresponding  to  this  rain  rate  distribution  from  (2)  Is 


N N 

A'  I A,4  I a(lU  dB  (9) 

1*1  1 n 1=1  1 

where  and  are  the  total  attenuation  and  rain  rate  for  the  1^ 
cell. 

The  general  model  of  (8)  and  (9)  has  unknowns  N and  { > . To 
simplify  the  solution  a ten-cell,  two-level  rain  rate  model  was  selected. 
The  unknowns  are  then  C,  the  fraction  of  the  path  over  which  the  rain  rate 
extends  and  which  is  a multiple  of  0.1L,  and  x,  the  parameter  which 
specifies  the  rain  rate  of  the  remaining  cells,  i.e. 


R 


x 


R-j  1 0 mm/hr 


R-j  >.10  mm/hr 


(10) 


The  attenuation  expression  (9)  then  reduces  to 


A = [C  afRj)  + (1-C)  a(Rx)]L  dB  . (11) 


Attenuation  and  rain  rate  data  pairs  for  the  experiments  Bl,  B8,  Cl,  C4, 
and  CIO  listed  In  Table  I were  used  to  determine  x and  C such  that  the 
predicted  attenuation  from  (10)  and  (11)  best  fit  the  measured  values 
over  the  range  of  rain  rates  for  the  experiments.  Specific  attenuation 
parameters  a^  and  bj  In  (1)  were  computed  for  the  elevation  angle, 
frequency,  polarization,  etc.  of  the  experiment.  [5]  The  best  fit  results 


- 17  - 


are 


IJ 


U 

II 

II 

li 

II 


‘ 


0 

a 

o 


D 

a 

i 

hj 


i 

$ 


i 


! 


-E 


x * -0.66  and  C = 0.2  . (12) 

This  result  was  also  obtained  earlier  [5,16,17]  by  examining  attenuation 
data  for  a single  experiment  (VPI&SU  COMSTAR  28.56  GHz  summer  1977  data) 
and  a single  rain  rate  ( R-|  =40  mm/hr).  We  also  note  that  a two-level 
rain  rate  distribution  for  terrestrial  paths  has  yielded  good  agreement 
to  attenuation  measurements.  [18] 

Theoretical  predictions  using  (10)  - (12)  and  the  values  of  L in 
Table  I are  compared  to  measured  data  in  Fig.  4 for  20  and  30  GHz  and  in 
Fig.  5 for  11  GHz.  Specific  attenuation  parameters  a^  and  b^  used  in  (11) 
were  developed  using  the  Rain  Propagation  Prediction  model  of  [5]  for  the 
polarization  state,  frequency,  and  elevation  angle  of  the  experiment. 

Also,  a mixture  of  40%  spherical  and  60%  oblate  rain  drops  with  a zero 
mean,  12°  standard  deviation  normal  distribution  of  canting  angles  was 
used  to  generate  a-j  and  b1 . Note  that  in  Fig.  4 agreement  is  very  good 
for  various  locations,  frequencies,  elevation  angles,  etc.  Predictions 
for  11.7  GHz  are,  however,  consistently  low  as  noted  by  the  data  from 
three  experiments  shown  in  Fig.  5.  This  discrepancy  is  possibly  due  to 
storms  which  pass  through  the  propagation  path  but  not  over  the  site  rain 
gauge.  This  effect  is  one  of  a lower  than  actual  rain  rate  for  a given  attenua- 
tion and  is  more  pronounced  at  11  GHz  than  at  20  or  30  GHz  because  of  the 
higher  rain  rates  required  for  noticeable  attenuation  at  11  GHz.  In  any 
event,  the  model  may  be  verified  at  11  GHz  by  examining  isolation  versus 
attenuation,  which  removes  the  point  rainfall  rate  dependence.  Experimental 
data  and  theoretical  predictions  [5]  based  on  the  two-level  rain  distri- 
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30.0  55°  5.493  Austin,  TX  [14]  30°  May  1975  ATS-6 

__  _______  COMSTAR 

28.56  44°  5.758  Blacksburg,  VA  [3]  37°  June  1978  D-2 
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Figure  4a.  Measured  data  from  the  Comsat  Labs. 

COMSTAR  D-l  19.04  GHz  experiment 
(Bl)  compared  to  predicted  results 
from  the  synthetic  storm  model 
(solid  curve). 
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Figure  4b.  Measured  data  from  the  NASA  Rosman 
ATS-6  20  GHz  experiment  (B2) 
compared  to  predicted  results  from 
the  synthetic  storm  model  (solid 
curve).  The  dotted  curve  is  a fit 
to  the  data  given  by  1.1547  Ro<6990. 
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Figure  4d.  Measured  data  from  the  Comsat  Labs. 

COMSTAR  D-l  28.56  GHz  experiment 
(Cl)  compared  to  predicted  results 
from  the  synthetic  storm  model 
(solid  curve). 
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Figure  4e.  Measured  data  from  the  NASA  Rosman 
ATS-6  30  GHz  experiment  (C3) 
compared  to  predicted  results  from 
the  synthetic  storm  model  (solid 
curve).  The  dotted  curve  is  a fit 
to  the  data  given  by  1.7748  R0,8081 
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Figure  4g.  Measured  data  from  the  VPI&SU  COMSTAR 
D-2  28.56  GHz  experiment  (CIO) 
compared  to  predicted  results  from 
the  svnthetic  storm  model  (solid 
curve) . 
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Figure  5a.  Measured  data  from  the  U.  of  Texas  CTS 
11.7  GHz  experiment  (A8)  compared  to 
predicted  results  from  the  synthetic 
storm  model  (solid  curve). 
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Figure  5c. 


Measured  data  from  the  VPI&SU  CTS 
11.7  GHz  experiment  (A13)  compared 
to  predicted  results  from  the 
synthetic  storm  model  (solid  curve), 
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butlon  model  are  shown  In  Fig.  6.  Note  the  good  agreement  for  all  three 
experiments. 

Recently  Hodge  [19]  has  proposed  a method  for  frequency  scaling  of 
attenuation.  Assuming  a Gaussian  spatial  distribution  of  rain  rate  the 
ratio  of  attenuations  AM  and  A'  at  frequencies  f"  and  f*  Is 


(13) 


where  R Is  the  maximum  of  the  rain  rate  distribution  and  a.j,  bj,  a!j, 
and  b^  are  the  power  law  fit  coefficients  for  specific  attenuations  as 
In  (1).  Using  the  specific  attenuations  for  30  GHz  and  20  GHz  at  an 
elevation  angle  of  44°  In  (13)  gives 


2.49  R"0*0528  . 


(14) 


The  synthetic  storm  model  results  yield 
jffjjj.  - 2.57  R-0-06'2  . 


(15) 


This  and  the  Hodge  model  of  (14)  are  plotted  in  Fig.  7.  Note  the 
excellent  agreement  between  the  two  the  etical  predictions.  Attenuation 
ratio  data  is  plotted  in  Fig.  8 for  the  28.56  and  19.04  GHz  data  from 
the  VPI4SU  COMSTAR  D-2  experiment  during  April  to  June,  1978,  as  a 
function  of  19  GHz  attenuation.  The  bars  indicate  the  standard  devia- 
tion of  the  data.  The  solid  curve  is  the  theoretical  prediction  using 
the  synthetic  storm  model.  These  attenuation  ratio  results  (both  theory  and 
data)  are  In  good  agreement  with  that  reported  for  ATS-6  30  and  20  GHz  measure- 
ments [14,20,21],  as  well  as  for  other  COMSTAR  28  and  19  GHz  measurements  [13]. 
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ISOLATION  VS.  ATTENUATION 
DATA  SET  A8 


Figure  6a. 


ATTENUATION  (DB) 


Isolation  versus  attenuation  data  from 
the  U.  of  Texas  CTS  11.7  GHz  experiment 
(A8)  compared  to  predicted  results  from 
the  synthetic  storm  model  (solid 
curve) . 
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Figure  6b.  Isolation  versus  attenuation  data  from 
the  VPI&SU  CTS  11.7  GHz  experiment 
compared  to  predicted  results  from  the 
synthetic  storm  model  (solid  curve). 
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Figure  6c.  Isolation  versus  attenuation  data  from  the 
Bell  Labs.  CTS  11.7  GHz  experiment  (A15) 
compared  to  predicted  results  from  the 
synthetic  storm  model  (solid  curve). 
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Figure  7. 


Ratio  of  rain  attenuation  at  30  GHz  to  that 
at  20  GHz  as  a function  of  rain  rate 
computed  from  the  synthetic  storm  model 
(solid  curve)  and  from  the  Hodge  model 
(dashed  curve).  The  elevation  angle  is  44° 


Figure  8.  Ratio  of  28.56  GHz  attenuation  to  that 
at  19.04  GHz  measured  at  VPIASU  from 
the  COMSTAR  D-2  satellite  during  April, 
May,  and  June  of  1978.  Error  bars 
Indicate  one  standard  deviation.  The 
solid  curve  Is  that  predicted  using  the 
synthetic  storm  model. 
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D.  ATTENUATION  PREDICTION  USING  RAINFALL  ACCUMULATION 
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The  least  expensive  method  for  predicting  rain  attenuation  statistics 
on  an  earth-space  link  would  involve  the  use  of  ground  rainfall  accumulation 
data  only.  This  approach  is  very  desirable  for  evaluating  earth  station 
locations,  because  widely  available  rainfall  accumulation  data  used  to 
predict  attenuation  statistics  would  preside  system  outage  estimates. 

The  general  solution  to  this  problem  is  in  its  infancy,  but  in  this  section 
the  synthetic  storm  model  will  provide  some  encouraging  results. 

The  Rice-Holmberg  rain  rate  model  [6]  provides  a method  of  predicting 
rain  rate  statistics  from  rain  accumulation  data.  It  yields  the  percent 
of  time  that  rain  rate  R is  exceeded  over  a time  period  P in  hours  during 
which  M millimeters  of  rain  fell  as  follows 

%T(R)  * RQ{ [0 . 03  6 e-°-03R  + 0.2(1  - e)[e"°*258R  + 1.86  e~L63R]}  (16) 

where  RQ  * 100  M/P  and  6 is  the  ratio  of  accumulation  of  "thunder  storm 
rain"  to  the  total  rain  accumulation.  The  attenuation  results  from  the 
synthetic  storm  model  can  be  fitted  to  a power  curve  giving  A = aR^  which 
in  turn  can  be  solved  for  R as 

R = expjV1  1n(§)]  . m/hr  (17) 

u 

Using  (17)  in  (16)  will  yield  the  percent  of  time  that  attenuation  A is 
exceeded . 

This  procedure  was  tested  with  statistical  rain  and  attenuation 
data  collected  at  the  VPI&SU  earth  station  for  the  calendar  year  1978. 

During  that  time  (P  = 8766  hr)  a total  rain  accumulation  of  M=434  im  was 
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measured.  The  fraction  of  thunderstorm  rain  6 was  chosen  to  be  0.2  from 
the  ratio  of  rain  accumulation  for  rates  above  30  irm/hr  to  the  total 
accumulation;  this  value  agrees  with  that  in  [6].  These  values  in  (16) 
yield  the  solid  curve  in  Fig.  9.  Also  shown  are  the  measured  rain  rate 
statistical  data.  Note  the  excellent  agreement  to  the  simple  Rice-Holmberg 
model . 

The  theoretical  predictions  of  attenuation  statistics  using  the 
predicted  rain  rate  statistics  of  Fig.  9 and  the  attenuation  relationships 
from  the  synthetic  storm  model  are  plotted  as  the  solid  curves  in  Fig.  10. 
The  measured  data  for  the  COMSTAR  (28.56  and  19.04  GHz)  and  CTS  (11.7  GHz) 
experiments  are  also  shown.  Agreement  is  very  good,  except  for  a slight 
shift  between  theory  and  experiment  for  CTS  which  was  noted  previously 
for  several  CTS  experiments. 
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Figure  10. 


ATTENUATION  (OB) 


Attenuation  exceedence  for  the  year  of  1978 
measured  at  VPI&SU  at  three  frequencies: 
28.56  GHz  (0),  19.04  GHz  (X),  and  11.7  GHz 
(+).  The  solid  curves  are  those  predicted 
from  the  Rice-Holmberg  rain  rate  of  Fig.  9 
and  the  attenuation-rain  rate  relationships 
from  the  synthetic  storm  model. 
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This  part  of  the  report  has  illustrated  that  the  effective  path  length 
model  for  predicting  total  attenuation  from  specific  attenuation  has 
limited  application.  The  synthetic  storm  model  was  proposed  as  an 
alternative.  The  piecewise  uniform  rain  rate  model  with  N cells  along 
the  rain  portion  of  the  coirmuni cation  path  appears  to  have  sufficient 
flexibility  for  a variety  of  modeling  applications.  Good  agreement  to 
experimental  results  was  obtained  using  only  two  rain  rate  levels  along 
the  path.  As  a larger  data  base  for  many  frequencies,  locations,  and 
elevation  angles  becomes  available,  a larger  portion  of  the  model 
flexibility  may  yield  further  improvements,  as  well  as  permitting  models 
for  particular  climatic  zones.  Experimental  confirmation  demonstrated 
the  tendency  of  the  model  to  be  independent  of  frequency,  location,  and 


elevation  angle. 
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It  was  also  shown  that  the  synthetic  storm  model  can  be  used  to 
scale  attenuation  from  one  frequency  to  another.  In  addition,  together 
with  the  Rice-Holmberg  rain  rate  model  the  synthetic  storm  model  can 
be  used  to  predict  attenuation  statistics  from  ground  rainfall  accumul l- 
tion  alone. 
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RECIPROCITY  AND  MILLIMETER 
WAVE  PROPAGATION  THROUGH  RAIN 


G.  INTRODUCTION 
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The  nonreciprocal  nature  of  millimeter  wave  propagation  through 
precipitation  has  received  little  attention  in  the  literature.  Nevertheless, 
in  some  applications  the  reciprocal  and  nonreciprocal  properties  of  a 
communication  system  with  precipitation  along  the  propagation  path  are 
important.  For  example,  the  downlink  signals  on  an  earth-space  conmuni- 
cation  link  could  be  used  in  a polarization  orthogonality  restoration 
system  [22] . The  downlink  effects  on  the  signals  must  be  the  same  as  on 
the  uplink  (i.e.  the  medium  must  be  reciprocal)  in  order  to  correctly 
apply  adaptive  cancellation  to  the  uplink.  Also  reciprocity  can  play  a 
role  in  millimeter  wave  experiments  involving  precipitation  media. 

In  this  part  of  the  report  we  discuss  those  characteristics  of  a 
propagation  medium  which  lead  to  nonreciprocal  behavior  and  the  signal 
parameters  which  they  affect.  A few  simple  examples  for  millimeter  wave 
signals  passing  through  rain  are  presented  together  with  numerical 
results. 
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H.  THEORY 

Consider  a millimeter  wave  signal  propagating  in  the  +z-direction  with 
a depolarizing  medium  (such  as  precipitation)  somewhere  along  the  path. 

The  output  electric  field  components  are  given  by 


■ [D] 


D D E 

xx  xy  x 


V Vi  LEy 


where  E ^ and  E 1 are  the  electric  field  components  associated  with  the 
x y 

incident  wave.  The  depolarization  matrix  [D]  represents  the  change  in 
the  incident  field  due  to  the  presence  of  the  depolarizing  medium;  it  is 
the  identity  matrix  for  an  entirely  free  space  path.  There  are  various 
models  available  [24,25,5]  for  determining  the  complex-valued  entries  of 
[D]  for  precipitation  media.  However,  for  this  investigation  into 
nonreciprocal  behavior,  the  general  properties  of  [D]  are  of  primary 
interest. 

If  the  received  signals  on  a dual  polarized  propagation  link  are  the 
same  upon  interchange  of  the  transmitter  and  receiver,  the  medium  is  said 
to  be  reciprocal.  The  well-known  reciprocity  theorem  states  that  if  a 
medium  is  linear  and  isotropic,  but  not  necessarily  homogeneous,  it  is 
reciprocal.  But  a precipitation  medium  such  as  rain  is  anisotropic  since 
it  contains  oblate  raindrops  with  scattering  properties  that  depend  upon 
the  Incident  polarization.  It  is  usually  assumed  that  oblate  raindrops 
are  either  equioriented  or  canted  symmetrically.  Then  reciprocity  still 
holds.  But  if  the  medium  is  also  inhomogeneous  along  the  propagation 
path  with  respect  to  canting  angle,  then  the  medium  will  be  nonreciprocal. 
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To  see  this,  consider  a piecewise  inhomogeneous  path  consisting  of  two 
cells,  each  of  which  Is  linear,  homogeneous,  and  symmetrically  anisotropic. 
The  homogeneous  condition  for  each  cell  permits  uniform  mixtures  of 
various  particle  sizes  and  shapes.  Let  the  depolarization  matrices  [C^] 
and  [Cg]  represent  the  cell  behavior  in  a manner  similar  to  (18).  Suppose 


an  incident  wave  with  components  E J and  eJ  enters  cell  1 first  and  then 

a y 


cell  2,  which  we  shall  refer  to  as  case  a.  Then 


E a E 1 

x - [C2][C,]  x 
a 2 1 E 1 

.y  J L_y 


The  same  input  wave  propagating  in  the  reverse  direction  (case  b)  produces 
the  following  output  field  components 


E b~j  [~E  1 

x ■ [CJ[CJ  x 
b ' 2 E1 

y _y 


The  overall  depolarization  matrices  for  the  two  cases  are 


[Da]  • CC23EC1]  and  [Db]  = [C^ECg] 


Performing  the  indicated  matrix  multiplications  in  (21)  and  comparing 
entires  in  the  matrices  [Da]  and  [Db]  reveals  that 


Da  = Db,Da  = Db, 
xx  xx  yy  yy 


Da  = Db,Da  = Db. 
xy  yx  * yx  xy 


(22a) 

(22b) 


This  result  Is  Independent  of  the  properties  of  the  cells,  except  that 
each  cell  must  have  a symmetric  matrix. 
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For  the  piecewise  Inhomogeneous  medium  to  be  reciprocal,  the  total 
path  effects  must  be  the  same  for  the  two  propagation  directions.  Thus 
[Da]  * [Db],  i.e.  [C-j]  and  [C2]  must  commute.  Equating  entries  from  the 
results  of  (21)  gives.  In  addition  to  (22), 

“xy8  * V3  »"•*  Dxyb  ' Dyxb  ' (23) 

Therefore,  the  depolarization  matrix  for  the  entire  medium  must  be 
symmetric.  This  is  a familiar  result  for  two -port  network  scattering 
matrices. 

If  the  total  path  depolarization  matrix  is  not  symnetrlc,  in  general 
signal  attenuation  and  crosspolarization  will  not  be  reciprocal.  Consider 
a precipitation  path  with  two  rain  cells  each  with  arbitrary  but  uniform 
rain  rate  and  canting  angle.  This  two  cell  medium  is  reciprocal,  i.e. 

(23)  is  satisfied,  if  and  only  if  [5] 


9„  - e-i  = n 


n = 0,1 


0 

a 


fi 

i 

i s 


where  and  e2  are  the  canting  angles  for  the  oblate  raindrops  in  the 
cells  1 and  2.  Thus,  nonreciprocity  arises  from  discontinuities  in  the 
canting  angles  of  the  rain  along  the  path,  in  particular  when  the  principal 
drop  axes  are  not  aligned  and  do  not  differ  by  90°.  Equation  (24)  also 
holds  for  two  cells  each  having  a symmetric  distribution  of  canting  angles 
where  and  e2  are  now  mean  canting  angles  for  the  distributions. 

The  attenuation  and  cross  polarization  introduced  by  the  medium  are 
of  interest  for  communication  system  calculations.  These  can  be  obtained 
by  computing  the  components  of  the  electric  field  intensity  exiting  the 
medium  which  are  co-  and  cross-polarized  to  the  incident  electric  field. 

For  arbitrary  incident  polarization,  attenuation  and  crosspolarization 
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ratio  (CPR)  are  not  reciprocal  if  the  raindrops  are  canted  differently 
along  the  propagation  path.  In  most  situations  the  effects  on  attenuation 
are  negligible  and,  in  fact,  it  can  be  shown  that  the  attenuation  is 
exactly  reciprocal  for  linear  input  polarization.  However,  in  certain 
situations  CPR  differences  for  opposite  propagation  directions  can  be 


signifi cant. 


I.  EXAMPLES 


As  a simple  example,  consider  a rain  medium  with  two  uniform  cells  of 
arbitrary  rain  rate  and  extent.  Let  cell  1 have  all  equioriented  oblate 
raindrops  aligned  with  their  major  axes  horizontal.  In  cell  2 all  drops 
are  aligned  but  differently  from  those  of  cell  1.  The  incident  polari- 
zation is  vertical  linear.  Then  the  difference  in  CPR  between  forward 
and  reverse  propagation  is 

ACPR(dB)  = CPRa(dB)  - CPRb(dB)  = A^fdB)  - Ayi (dB)  (25) 

where  A^(dB)  and  Ay^(dB)  are  attenuations  in  dB  of  cell  1 for  linear 
polarizations  along  the  major  and  minor  drop  axes,  respectively.  [5] 

Note  that  this  result  is  independent  of  the  canting  angle,  length  and  rain 
rate  of  cell  2.  For  the  case  of  a 2 km  length  of  cell  1 along  a path  with 
a 40°  elevation  angle,  the  CPR  difference  of  (25)  is  plotted  in  Fig.  11 
for  three  frequencies  using  available  differential  attenuation  values  [23]. 
Note  that  the  CPR  difference  increases  with  increasing  rain  rate. 

Another  example  illustrates  the  effects  of  different  canting  angles 
in  the  two  cell  rain.  Again  the  incident  polarization  is  vertical  linear, 
but  the  canting  angles  (of  drop  major  axes  relative  to  horizontal)  of 
drops  in  cells  1 and  2 are  and  62*  For  this  case,  and  for  an  identical 
rain  rate  and  path  length  in  both  cells, 

II  + g cot  e,  tan  e9l 

ACPR(dB)  = 20  log  '■ (26) 

|g  + cot  0-j  tan  e2| 

where  g a expta^  - ay  + j(BH  - ey) ]L  in  which  a^  - ay  and  bh  - By  are  the 
differential  attenuation  and  phase  constants  and  L is  path  length  of  each 


R^ (mm/hr) 


5i 


Figure  11 


. Difference  in  CPR  for  opposite  directions  of  propagation 
through  two  rain  cells  as  a function  of  rain  rate  in 
cell  1,  which  has  a canting  angle  of  0°  and  a length  of 
2 km.  Cell  2 has  arbitrary,  but  nonzero,  rain  rate, 
canting  angle,  and  length.  The  Incident  polarization 
is  vertical  linear  and  the  path  elevation  angle  is  40°. 


cell.  [5]  The  CPR  for  forward  and  reverse  propagation  are  shown  in  Fig.  12 
for  30  GHz,  20°  elevation  angle,  cell  lengths  of  5 km,  a rain  rate  of 
25  mm/hr  in  both  cells,  0-j  = 15°,  and  e2  variable.  Note  that  ACPR(dB)  = 0 
for  e2  - 15°  when  all  drops  are  aligned  when  n = 0 in  (24),  and  for  e2  = 105° 
when  n = 1 in  (24).  The  CPR  difference  is  also  zero  for  e2  = 165°  when 
e2  * it  - and  the  drops  are  symmetrically  canted  having  a cancelling 
effect  in  this  case. 
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Figure  12.  Crosspolarization  ratio,  CPRa  and  CPRb,  for  forward  and  reverse 
propagation  and  their  difference,  |aCPR|,  for  two  rain  cells 
each  of  length  5 km,  as  a function  of  the  canting  angle  of 
cell  1.  The  canting  angle  of  cell  2 Is  15°.  The  rain  rate 
for  both  cells  is  25  m/hr.  The  frequency  is  30  GHz,  the 
elevation  angle  is  20°  and  the  incident  polarization  is 
vertical  linear. 
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